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Abstract 



, It is shown that only a narrow class of inflationary models can possibly 

agree with the available observational data on the anisotropy of the cosmic 



microwave background radiation (CMBR). These models may be governed 
by "matter" with the effective equation of state —1.2 < p/e < —0.6 which 
includes the De-Sitter case p/e = —1. 
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The recent discovery of the angular variations in the CMBR [1] has strongly sharpened 
the issue of the origin and nature of the long-wavelength cosmological perturbations. The 
consequences of inflationary models are under active investigation (see, for instance, a recent 
paper [2] and references therein). It seems to me that only a narrow class of the inflationary 
models, discussed below, can possibly avoid theoretical and observational inconsistencies. 

I am considering here relic gravitational waves. The variable gravitational field of all 
cosmological models (unless the cosmological scale factor a(rj) is such that a" = 0), and 
inflationary models as a particular case, inevitably generate gravitational waves [3]. (The 
graviton creation in FRW fields was denied in the past and anisotropic models were claimed 
to be the only way to get a nonzero result, but now they are not considered anymore as a 
necesary condition.) The generating mechanism is quantum- mechanical in its nature, and 
the generated perturbations are always placed in squeezed vacuum quantum states (see [4] 
and references therein). This means that different modes of the created field are not totally 
independent, as is often assumed in the literature on inflation, but, on the contrary, some 
of them are highly correlated which leads to the picture of standing waves and modulated 
spectra. The generated gravitational waves inescapably produce the angular anisotropy 
in the CMBR. The angular correlation function for 5T/T variations caused by squeezed 
gravitational waves has been derived recently [4]. We will use it here in our analysis. 

The dimensionless gravity-wave field, with all the normalization factors taken into ac- 
count, can be written as 

/ , 1 roo 2 i r . , 

M^x) =4^^--— _ d 3 n£^.(n)-= [c£fo) e *« + c^e— ] (1) 



a( V ) (2tt) 3 / 2 J -co £T« v 'y/frt 
where the scale factor a(i]), ds 2 = a 2 (r))(dr) 2 — dx 2 — dy 2 — dz 2 ), has the dimension of 
length and l pi = (Gh/c 3 ) 1 ^ 2 is the Planck's length, all other quantities are dimensionless. 
The two "transverse-traceless" polarization tensors p^-(n) (s = 1, 2) satisfy the conditions 
pfj (n)p s '' lJ (n) = 25 ss ', pfj(-n) = pfj(n). The time-dependent annihilation and creation 
operators c^(r]), c„ (r?) can be written (for each s) as 

c n (v) = u n (rt)c n (0) + v n (rj)cl n (0), 4 (77) = <fa)4(0) + v* n ( V )c_ a (0), (2) 



where c n (0), (4(0), are the initial values of the operators taken at some i] = i] long before 
the interaction became effective. The complex functions u n (rf), v n (i]) satisfy the equations 

iu' n = nu n + i(a! / a)v*, iv' n = nv n + i(a' / a)u* n (3) 

where ' = d/drj, \u n \ 2 — \v n \ 2 = 1 and u n (0) = 1, v n (0) = 0. It follows from these equations 
that the function fi n (r]) = u n (r)) + obeys the equation ji" n + (n 2 — a" /a)ji n = which 

is precisely the equation for classical complex /i-amplitude [3] of the gravity-wave field. 

The (Bogoliubov) transformation (2) can be written in a form involving the two-mode 
squeeze operator which demonstrates the inevitable appearance of squeezing in this kind of 
problems. In the Schrodinger picture, the initial vacuum state |0), c n (0)|0) = 0, evolves into 
a strongly squeezed vacuum state. 

The angular correlation function for S/T (see Eqs. (12), (13) in Ref. [4]) can be rearranged 
by using the "summation theorem" [5] and the formulas relating the Gegenbauer polynomials 
to the associated Legendre polynomials [6] and cast into an elegant exact form which directly 
involves the Legendre polynomials P[(cos5): 

AT AT 00 
<0|^r(e*) — (e*)|0> = K^oxS) (4) 
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where K\ = (21 + 1)1(1 + 1)[Z(Z + 1) - 2\F h 
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and 



fn(VR ~W) = -^=(fJL n /a)' . 

v2n 



The two unit vectors e\, e\ point out in the directions of observation and 8 is the angle 
between them. The photons of CMBR are assumed to be emitted at 77 = r] E and received 
at (present) time 77 = r] R ; w = r/ R — 77, w\ = rjR — r/ E (the Sachs- Wolfe effect [7]). Note that 
the correlation function, with no additional assumptions made whatsoever, is rotationally 
symmetric and its multipole expansion begins from 1 = 2. 



The derived formula is universal and can be used with arbitrary a(r)). We will apply this 
formula to simple models consisting of three consecutive stages of expansion: inflationary 
(i-stage), radiation-dominated (e-stage) and matter-dominated (m-stage) [8]. 

The scale factor of the model can be written at the three stages as follows, 
i-stage: 

a (v) = ^oM 1+/3 , V < Tj! , r?i < 0, 

e-stage: 

a(rj) = l a e (rj - rj e ), Vi < V < V2 

where a e = -(1 + / 9)|t7 1 | /3 , r] e = j^Vi, 
m-stage: 

a(rj) = l a m (rj - r? m ) 2 , r] 2 < rj, 

where a m = [a e /4(?72 — r) e )], r\ m = —r] 2 + 2r] e . The functions a(i]), a'(r]) are continuous at 
rj = r]i and r] = 772. All expanding models with 1 + f3 < (rj must be negative if 1 + f3 < 0) 
are inflationary in the sense that the length scale equal to the Hubble radius at some early 
time of expansion can grow at all three stages up to, at least, the size of the present day 
Hubble radius l H = a 2 /a', rj = rj R . The case (3 = — 2 corresponds to the De-Sitter expansion, 
the cases (3 < — 2 correspond to the so called power-law inflation (a(t) ~ t m , m > 1) and 
the cases —2 < (3 < —1 (apparently, not having been analyzed before) correspond to the 
law of expansion a(t) ~ \t\ m , m < —1, t < 0. The i-stage is governed by "matter" with the 
effective equation of state p = q((3)e, where q(f3) = (1 — /3)/3(l + (3) and q(f3) varies from 
— 1/3 to —00 for —00 < (3 < —1. Expansion is accompanied by the growth of energy density 
and curvature if— 2 < (3 < — 1. 

In realistic cosmological models a(r)E) /a(r)R) « 1CT 3 , a(rj 2 ) / a(r] R ) as 1CT 4 . Also, one 
has 3.1CT 32 < aijji) / a{rj R ) < 3.1CT 12 if one wants to commence the e-stage at densities not 
lower than the nuclear and not higher than the Planckian, or, in other words, if one wants 
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the Hubble radius k at the end of inflation, U = — l a (l + f3) 1 |?7i| 2+/3 , to be in the interval 
1 < k/l pl < 10 40 . 

To define the numerical values of i] it is convenient to choose t)r — r\ m = 1. Then, 
« 5- 2 //3|i + /3|-V/3 (Z /Z H )-V/» and fe// ) 2 » (25^// H ) 2+/3 |1 + 0| 2(1+/3) ( The 
wavelength A = 2na/n equal to Z# has the wave number n H = 4ir, and the wavelength equal 
to the Hubble radius at rj = r\ 2 corresponds to n m ~ 4ir • 10 2 . The minimally sufficient 
inflation should begin not later than at r] b = (1 + /3)/2, its variable gravitational field 
generates waves with wavelengths up to In- Inflation that started earlier inevitably generates 
the longer waves also. 

The gravity- wave eguation using the scale factor of the form a{rf) = a r] 1+f3 has been 
solved and the relation between the initial and final amplitudes has been derived earlier [3] . 
The waves start oscillating with the amplitude B{n) which is related to the initial amplitude 
A{n) by B 2 (n) ~ A 2 (n)(nr)i) 2 ( 1+ ^ (one can use formula (5b) from the second Ref. [3] in 
which the interpretation of the participating amplitudes should be reversed because, for 
1 + (5 < 0, the condition (ni]) 2 ^> 1 is satisfied initially and gets violated later on). Since 
the initial (vacuum) spectrum goes as A(n) ~ n this leads to B(n) ~ n /3+2 and B(n) ~ n° 
for p — — 2, that is, in case (3 — — 2, all waves start oscillating with the same amplitude (the 
Harrison-Zeldovich spectrum). In case of (3 < — 2 the spectrum gets "tilted" by increasing 
the relative contribution of longer waves, and in case of (3 > — 2 — by increasing the relative 
contribution of shorter waves. 

The exact solution to Eq. (3) for the complex function (J, n (r)) satisfying the required 
initial data and continuous with its first time-derivative at the joining points rji, rj 2 has the 
following form, 
i-stage: 

Unfa) = (nv) 1/2 [ A iJp+i/2(nri) + A 2 J-^ +1 / 2 )(nrj)] 

(for a technical simplification we work solely with the Bessel functions and exclude the 
half-integer /3's but the final result will be free of this limitation) where 
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Ai = l -r\l^e^ 0+ ^ 2 \ A 2 = iA 1 e-™ (3 and x = n Vo , n\r) \ > 2tt|1 + /3| 

COS D7T V 2 



COS /3-7T 

e-stage: 



my 



where y = n(r) — r) e ), 
m-stage: 

Mnfa) = ^(Cl + C 2 J_ 3 /2(^)) , (6) 

)' = (Ci - C 2 J_ 5/2 (z)) (7) 

where z = n(?7 — r] m ). Note that J-z/ 2 (z) represents the so-called decaying solution which is 
necessarily present. 

The coefficients Bi, B 2 , Ci, C 2 , are determined by the continuos joining of the solutions. 
In particular, 

-d = B x {a 2 + (5*) + B 2 (a* 2 + fo) , -iC 2 = B x {a 2 - (%) - B 2 (a* 2 - fa) 

where a 2 = e iy2 (8y 2 — 1 + iAy 2 )/8y 2 , (3 2 = —e l3y2 /8y 2 , y 2 = n(?? 2 — r] e ). We are only 
interested in modes that have interacted with the barrier U{rj) = a" /a and have been 
generated quantum-mechanically. Their wave numbers obey the condition n|?7i| <C 2tt\1 + (3\ 
and for them 

B x » -B 2 » i e "°(/5 + 1)^( / 9)H 1 )' 3 = 5, 

where 

= y|e i7r/3/2 [cos/37r2 /3+1/2 r(/3 + 3/2)] _1 , = 1 for (3 = -2 . 

The values of Ci, C 2 depend on whether ?/ 2 > 1 or ?/ 2 < 1. For relatively short waves, 
n ^> n m , one has approximately C\ ~ — 2iBsiny 2 , C 2 ~ 2iBcosy 2 . For longer waves, 
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n <C n m , one has C\ « —^-By^ 1 , C 2 ~ — C 2 <C Ci. The additional large factor y^ 1 
in Ci reflects additional amplification of waves at m-stage. 

For a qualitative description of amplitudes and spectral slopes we introduce "charac- 
teristic" spectral components of the field: h(n) = l p in\fi n \/a, and its first time-derivative: 
h'(n) = lpi(\/jLn\/o)'. We have h(n) ~ ^n' 34 " 1 sm[n(r] — rj e )], h'(n) ~ l ^-n^ +1 cos[n(r] — rj e )] for 
n 3> n rn , and h(n) ~ l -^n^ cos[n(r] — r) m )], h'(n) ~ ^n 13 sm{n(r] — r) m )\ for n H ^> n ^> n m . 
For n % one can use the approximation z 1 in Eqs. (6), (7) and obtain 
h(n) w ^|^(/3)|^+ 2 , /i'(n) w ^|^(/3)|^+ 2 . 

We will now start deriving restrictions on inflationary models. Astrophysically interesting 
and consistent values of h require h(nn) to be not much larger or much smaller than 10~ 4 . 
The mean square value of the field diverges in the limit of small n {i.e. wavelengths much 
longer than l H ) for all (3 < -2 : (h 2 ) w (l pl /l ) 2 \ilj({3)\ 2 J ^ n ? +2 . This does not allow the 
duration of inflation to be arbitrarily long. In case of (5 = — 2 the divergence is logarithmic 
and the restriction on duration is very mild but it becomes increasingly severer for (5 < — 2. 
For (5 ~ —7, even the minimally sufficient duration of inflation does not help as the h{n H ) 
strongly exceeds 1CT 4 even if the largest allowed l i: k w 10 40 l ph is chosen. On the other hand, 
for (3 > — 2 there is no problem with the long-wavelength divergence but the amplitudes very 
quickly become too small. For (3 ~ —1.8, the h(nn) becomes smaller than the required level 
even if the smallest allowed /«, U ~ l p i, is chosen. (One may argue that the values U <C l p i 
do not necesserily invalidate the analysis and can also be allowed, although they imply the 
over-Planckian densities, since the wavelengths of our interest are much longer than l p i all 
the way up from the beginning of inflation. However, we will confine in this paper to the 
requirement U > l p \ which leads to (3 < —1.8.) 

Additional restrictions come from the ST/T considerations. We will be interested in the 
lower multipoles /, / < 30, to which the waves n ^> n m give a negligibly small contribution. 
Neglecting the second term in Eq. (7) one can write Eq. (5) as 



2 

,J i J l+1 / 2 (x) J 5 /2{n - x) 

x 5/2 ( n _ 2^3/2 



dn 



where wi ~ 1 — 3.10 2 . In the limit of long waves one can use the small argument approxi- 
mation for the Bessel functions and write 

where (f)(1) = [15^2 l T(l + l/2)l(l — l)] -1 . The quadrupole component K 2 is divergent in the 
limit of small n for all (3 < —4 so these values of (3 should be excluded unless the duration 
of inflation is precisely tuned. (Interestingly enough, the minimally sufficient inflation may 
be a likely outcome of certain quantum-cosmological models [9].) 

The multipole distributions K\ for some models have been numerically computed by A. 
Wiseman. For each /, the main contribution to Fi comes, as one might expect, from the 
interval of n between n ~ I and n ~ 21. In particular, K 2 is dominated by waves longer 
than Ijj [10]. The absolute values of K\ depend on (3 and l Q . For instance, the quadrupole 
component 10 2 (l o /l p i) 2 K 2 has the following numerical values: 9.7, 7.3, 5.6, 7.4, if one takes 
(3 = —1.8, —2.0, —2.4, —3.0, correspondingly. The multipole distributions normalized to 
K 2 , that is the functions Ki/K 2 , are independent of l Q . They are shown in Fig. 1. For 
completeness we have included also the extreme cases (3 = —1.6, —1.1, though they would 
imply strongly over-Planckian densities at the end of inflation, as was explained above. It 
is worth noting that the quadrupole component exceeds other multipoles and K 3 /K 2 < 0.6 
for all models with (3 < —1.8. The results for the (3 = — 2 case are in qualitative agreement 
with those in Ref. 11. It remains to be seen which of these distributions can survive after 
comparison with the detailed COBE-type observations. One should bear in mind, of course, 
that predictions for K\ are statistical and should be augmented with variancies based on the 
statistics of squeezed quantum states. 

In conclusion, the inflationary models governed by "matter" with the effective equation 
of state —1.2 < p/e < —0.6 seemingly avoid theoretical difficulties and some of them can 
possibly withstand comparison with the observations. (The models include the De-Sitter 
case if the duration of expansion was not excessively long.) Taking into account the density 
perturbations can only make this interval narrower. 
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FIGURES 



FIG. 1. The normalized multipole distributions. 
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